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We have earlier constructed a generalized entropy concept to show the direction of time in an 

evolution following from a Markov generator. In such a dynamical system, the entity found changes 

in a monotonic way starting from any initial state of the system. In this paper, we generalize 

the treatment to the case when population is pumped into the system from levels not explicitly 

considered. These populations then pass through the coupled levels and exit by decay to levels 

outside the system. We derive the form of the equation of motion and relate it to our earlier 

-—.. treatments. It turns out that the formalism can be generalized to the new situation. Its physically 

I relevant features are demonstrated, and the behaviour obtained is illustrated by numerical treatment 

/--s of the standard two-level system with pumping and relaxation included. 
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I. INTRODUCTION 



All irreversible evolution models should carry some signature of the direction of time. In the theory of dynamical 

systems [T] such an entity is called a Lyapunov function. This defines uniquely the forward direction of time. In 

r^ a physical context, the function may be considered to be a generalization of the entropy concept. In the case of 



approach to thermal equilibrium, the description should consistently reproduce what we take to be known from 

j^ thermodynamics. The procedure for finding this description will, of course, depend on the model we offer for the time 

^ evolution. The general concept of the direction of physical time is discussed rather comprehensively in [5] . 

O^ The modern theory of kinetic equations is compiled in the text [S]. The quantum counterparts are presented in 

[4] . Recent progress in quantum dynamics has created a separate research field named open systems dynamics [5] ; its 

^_^ more mathematical aspects are presented in [B]. 

^ The contemporary focus on information measures has led to the introduction of a multitude of generalized entropies 

^>0 [7] and [H]. The various expressions have different properties and different uses. In our earlier work, we have attempted 

^N to describe the evolution of a general dynamical system by constructing a state functional displaying the direction of 
^^ time by monotonically changing. This may then be regarded as a generalized entropy or, alternatively, as a Lyapunov 

function. 
'/^ In general, the time evolution equation of a physical system is of the first order and generated by a linear operator 

on the object defining the state of the system. For convenience, the generator is here taken to be time independent, 
which restricts the applicability of our approach. The time step is taken to depend only on the state at the present 
7^ time. This is usually called a Markovian evolution in physics. Within these restrictions we have been able to develop 
^ a rather detailed theory for the general case. 

In [9] the solution was presented for a restricted class of evolution equations. In the case of a driven two-level 
system with decay, the solution was presented in [ID!. As the cases discussed that far assumed no degeneracy, this 
d was treated in |llj . The case of a thermal reservoir was introduced in [12J . 

The method introduced by us displays a monotonically evolving quantity. This serves to define a forward direction 
of time. In our paper ^T3] , we derive the proper time-inversed behavior of irreversible time evolution. The related 
minimum property of entropy production is discussed in |14j . 

In this paper, we want to apply the theory to a case not usually addressed by the formal evolution theory. We 
consider the case where the physics takes place in a state space pumped by externally induced means. The population 
in the space considered is thus increasing with time, but this is counteracted by decay channels out of it. This 
system is not of the canonical Lindblad form [15], but it retains the physical interpretation by letting the density 
matrix elements denote the density of active particles instead of the customary probabilities. In laser physics and 
nonlinear spectroscopy this system has frequently been taken to model the physical situation investigated |16j . As 
most experimental laser physics is still described by equations of this class, we regard it to be important to relate 
such equations to the more recent approach developed by us. 

In Sec. 2, we introduce the concept of a pumped system and its quantum mechanical description by a linear master 
equation. In Sec. 3 we discuss the physical situation as it differs from a genuine Lindblad case. Finally we apply the 



theory introduced earlier to this case, and extract its characterization of the direction of time. The theory is then, 
in Sec. 4, illustrated by applying it to the case of a pumped decaying two-level system. This model has played a 
central role in both laser spectroscopy and quantum optics [T7]. The formalism, however, also applies to Markovian 
rate equations [T^. Finally the conclusions from the treatment are discussed in Sec. 5. 

II. THE PUMPED SYSTEM 
A. The time evolution equation 

The terminology of quantal state evolution assigns the term open system to any quantum system in contact with an 
environment. Ordinarily the system is probabilistically closed, i.e. we assume that the density matrix of the system 
retains the probability interpretation and is normalized. Consequently all probability exchange between the states, 
reversible as well as dissipative, occurs between levels in the system. The general memory-less expression for such 
evolution is given by the Lindblad form of the generating equation. 

There is, however, a class of time evolution equations that do not conform to this model. There is a flow of 
probability through the system; the states are pumped by some technical method and the influx of probability is 
compensated by states decaying to unobserved lower levels. Such models are termed ■pumped system,s here. Naturally 
the state space could be enlargened so that all pumping and decay mechanisms are included in the dynamics of the 
evolution. It is, of course, always possible to consider a space large enough to eliminate the need for a phenomenological 
master equation, but it is not always an expedient or even a practical possibility. 

The time evolution equation will now hold several different pieces: 

1. The decays inside and out of the system are described by a relaxation operator TZ, which we assume to be linear 
and memoryless. The time derivative of the density matrix in the absence of pumping and with H = 0, where 
H is the Hamilton operator, would be 

dtp = np. (1) 

2. The pumping is described by a random process which at time to introduces a system in the state 

n 

where {| (pn)} is an arbitrary basis on the state space of the system under consideration. At a later time this 
will, in the absence of relaxation, evolve to the state 

\^Pit))=J2r^n{t-to)\^n)- (3) 

n 

We now form an ensemble accumulated by all systems introduced at times ip l£ t- This is described by the state 

t 

Pnm{t) = / Cn{t-to)c*^{t~to)P4,{to)dto, (4) 

— CO 

where P,/i(io) is the probability of introducing a system in state ^p into the ensemble at time ig- This is assumed 
to be normalized over the states introduced: 

^P4t) = l. (5) 

The contribution from Q to the time evolution is given by two terms, 

t 

dtPnm{t) = Cn{Q)c*^{Q)P-^{t)+ J Ot [Cn {t - to) C*^ {t - to)] P^{to)dto. (6) 

— C30 

We now introduce the density operator averaged over the pumped states in ^ by setting 

y^ Pnm (t) -> Pnrn (t) ■ (7) 



We may then write the averaged time evohition as 

dtPnrn{t) = A„™(t) + Cnm{t), (8) 

where the inhomogeneous term is given by 

A„™(<) = ^c„(0)c:;,(0)P^(i). (9) 

Because the Hamiltonian is taken to be the same for all systems pumped into the ensemble, we let the state 
coefRcients satisfy the ordinary Schrodinger equation in the form 

dtCn {t - tf)) = -i^ HnkCk {t ~ to) , (10) 

where H is the generator of the unitary time evolution inside the system. Using this in (|6| we obtain 

C„™(t) = -»([iJ,p(i)])„„, (11) 

where we have introduced the ensemble density matrix from Hh. 



Now, collecting terms from (fTl) and (Is]), and using (111, we find the time evolution equation for the pumped ensemble 
to be 

dtp^A{t)~i[H,p{t)]+np{t). (12) 

This is the most general linear and memoryless evolution equation we encounter. The Lindblad form falls under this 
more general class if we set A{t) = 0. If this is not the case, then the relaxation term TZp cannot be of the Lindblad 
form, which on its own would ensure conservation of the probability. As the term A does not do this, its influence 
must be compensated by TZp. 

Not being of the Lindblad form, Eq. (12) does not manifestly generate completly positive time evolution. However, 
properly used, it does not violate the physical interpretation of the density matrix. The pumping in Eq. (9) is of the 
form of a density matrix and thus cannot cause any troubles. The damping operator (1) has to be introduced such 
that no abnormal results emerge. This implies inequalities between the decay constants, but these are no different 
from the corresponding relations encountered in Markovian rate processes, Ref. [18]. The Hamiltonian evolution, 
naturally, induces only physically acceptable changes. 

B. Properties of pumped evolution 



The equation ( 12 ) may be written in a notation generalizing the Lindblad case as 

dtp = A + Cp. (13) 

In agreement with the concept of pumping a state, it follows that 

A„„ = ^ c„ (0) < (0) P^ (t) > 0; (14) 

this term can only add to the population of the state. The pumping mechanism is usually incoherent, and A is not 
expected to have any non-diagonal contributions, because of the ensemble average. 
Conservation of probability requires that the steady state, if it exists, satifies 

Tr (Upo) = -TrA. (15) 

The operator TZ is thus essentially a negative operator. 
If no pumping occurs, A = 0, there are two possibilities: 

• After an infinite time, all population decays out of the system and the only steady state is the one with all 
populations vanishing, pq{oo) — 0. In this situation, the relation 

Cp^O (16) 

has got no nonvanishing solutions; all eigenvalues of C are nonzero. 



• There is a part of C which allows for a steady state with nonvanishing density matrix 

/:oP° = 0; pV 0. (17) 

This concerns a part of the state space which has no decay channel out of it, and population ending up here will 
remain trapped. Thus no part of the pumped population can go into this subspace, because then probability 
would accumulate here without limit. It thus describes a part of the state space which is totally decoupled 
from the rest of the system: No coherent transfer nor any irreversible decay can couple this to the rest. It is 
consequently a system which can be treated by the ordinary theory; all its population derives from its initial 
state and no pumping will affect it. In the subsequent discussion we omit this possibility. 



The steady state solution of ( 13 ) becomes 

Po = -C-^P^. (18) 



According to what is said above, the operator C has got no zero eigenvalues, and the expression ( 18 ) will exist in all 
physically acceptable situations. 

The ensemble density matrix elements p„„ > all denote the population on level n but they no longer necessarily 
relate to any probability; in particular it need not hold that 

Trpo^Trp{t). (19) 

Note that the special initial condition, p(to) = 0, is perfectly acceptable and may still imply po = p(oo) 7^ 0. The 
off-diagonal elements of p retain their role as determining the multipole moments of the system. 



With this notation we may write ( 13 1 as 

dtp ^C{p~po); (20) 

with the notation 

Sp = p-po. (21) 

This allows us to write 

dt6p = C5p, (22) 

which relates the present case to our earlier theory. 

III. PHYSICS OF TIME EVOLUTION 
A. Formal properties 

The objects p representing the state of the system obey the general dynamic time evolution equation of the form 

dtp = dtSp = A + Cp = CSp. (23) 

They can be taken to belong to a linear manifold of elements denoted by | p)). This space can be equipped with a 
natural inner product by writing 

{{p,\p2))^Tr{p\p2). (24) 

The time evolution operator is not assumed to be Hermitian with respect to the inner product defined. Consequently 
we need to introduce right eigenvectors 

£ I X,)) - A, I X,)) (25) 

and left eigenvectors 

cUy.)) = K\y.))- (26) 



From these relations follows that 

((2/m I ^ I ^^)) ^ K{{yt, I 2^1.)) 

(27) 
{{C^y^, I x^)) = Xf^iiVf, I x^)). 

Thus if Xi, 7^ A^ the states are mutually orthogonal and may be normalized against each other 

{{y, \ X,,)) = 6,^. (28) 

We assume in the following that both the right and left eigenstates form complete basis sets. This is, in particular, 
true if their corresponding eigenvalues are nondegenerate, but in general they are complex. 
With these definitions, we can present spectral representations for the operators 

(29) 

When both sets of eigenstates are complete, the identity operator can be written as 

3 = ^|2:,))((y. h^|y,))((x, |. (30) 

I' V 

The inner product between states belonging to the right eigenvalues have no simple relations. In order to obtain a 
simple situation, we introduce the mapping 



by setting 



The inverse mapping is seen to be 



f):{|x.))}^{|y,))} (31) 



f^ = El2^-))((y-l- (32) 



f7-i=E|x,))((x. |. (33) 



These are positive Hermitian operators, and it is possible to define a new metric based on the bilinear form defined as 

Mo[pi,P2] = ((Pi I ^ I P2)) = Tr{p\^p^). (34) 

This has got all the properties of an inner product, and it thus defines a topology in the space of all quantum states. 
Within this product we have 

Mn[|x,)),|x^))] = ((a;, |y,,))=<5,^. (35) 

A similar construction is possible on the states {| yv))^ by the use of fi^^.The operator fJ may be considered as a 
metric operator on the manifold of physical states. 
By direct calculation we find that 

nm-^ =Y,\y,))K{{x,\=C*K (36) 

Here as in the following, we assume that complex conjugation denoted by * affects only c-numbers, not states. 



From ( 36 ) , we derive the relations 



(37) 



forming the basis for the considerations below. 



If C has real eigenvalues, then C* = C and the operator is Hermitian with respect to the new inner product Mn 
defined above. If the eigenvalues occur in complex conjugate pairs, i.e. for each | x^)) there exists an | x*)) such that 



C\x*;)) = \l\x*^))^ (38) 

we find from 

Ap ((yp I x*^)) = {{y,, I C I x;)) = X; {{y, \ x*^)), (39) 

that 

{{y^ I xp) EE {{x^ I n I x;)) = 0; (40) 

these states are thus orthogonal in the Mq metric. We expect the eigenvalues of physical evolution operators to 
correspond to damped oscillations at frequencies given by the energy differences in the system. We call this case 
physical evolution. 



The time dependent state that is a solution of the evolution equation ( 23 ) is given by 

\5p(t)))=Y.'^ue^^{Kt)\x,)), (41) 

where 

r, = ((y, |<5p(0))). (42) 

From the physical meaning of the time evolution operator in a pumped system, we derive the following conclusions: 

• There exists no steady state | xq)) with eigenvalue Aq — 0. Asymptotically, the time evolution will lead to 
p{oo) = po. 

• Any solution must approach this in a smooth manner, thus we require Re A^ < for all v. 

B. The direction of time 

As the system undergoes irreversible time evolution, it singles out one direction of time defining the forward progress. 
In classical dynamics, we can achieve this by introducing a monotonically changing variational function; in system 
theory this is called a Lyapunov function. 

We find easily that the definition 

Mn{5p, 5p) = Tr {Sp'^ndp) = {{5p\n\5p)) (43) 

provides a function changing monotonically with time. We have namely 

dtM^ [{p - Po) , {P - Po)] - Tr (dtSp^nSp) + Tr (Sp^ndtSp) 

= Tr (Sp^C^nSp) + Tr {5p^nC5p) 

(44) 
= Tr{5p'^n{C*+C)5p) 

= 2Rc Tr {Sp^QCSp) . 



where we have used ( 37 1 . Using the expansion 

IMi)))=^^.exp(A,i)|a:^.)) 

we obtain 

dtMniSp, Sp)^Y.\ ^- I' (^*- + ^-) '^^P [(^- + ^-) ^] ^ 0' (45) 



because all eigenvalues have negative real parts. The rate (45) decreases with time and when steady state is reached, 



the development ceases. This is thus a proper variational operator which determines the direction of time towards 
the eventual steady state. In particular we note that the pumped system, in general, lacks a concept of thermal 
equilibrium. 



The expression ( 43 ) represents the natural extension of the purity to the case of irreversible time evolution. If we 



want an entity extensive in the combination of uncoupled subsystems, we may introduce the entropy-like expression 



5o = ±logMa; 
the choice of sign depends on the interpretation of the quantity. 



(46) 



IV. APPLICATION: PUMPED TWO-LEVEL SYSTEM 



We apply the considerations above to the generic driven two-level system with pumping and decay out of the system. 
The time evolution equation is now 
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Here lo is the detuning of the driving field, V is the Rabi type coupling between the levels, and F and 7 are the 
population and coherence decay rates. These define the time constants Ti and T2 respectively. The coherence decay 
satisfies the constraint 



7>^(ri + r2). 

This follows if we set y = A = and impose the condition 

Pll{t)p22{t) > Pl2{t)p2l{t). 



(48) 



(49) 



A. Analytic results 



It is straightforward to obtain the complex eigenvalues of the matrix in ( 47 ) and verify that all real parts correspond 



to decay. The steady state of (471 is also easily computed, but the ensuing expression is too unwieldy to allow an 



intuitive interpretation. One quantity of interest is the asymptotic population difference 

'A2 Ai\ (^ 7?2X/2 



^22(00) -pii(oo) 



Ti 



1 



a;2 +72 + j^2y2 



(50) 



The deviations from the pumping difference without coupling is in the form of a power-broadened Lorentzian, where 
the rate of incoherence is measured by the dimensionless parameter 



^2 ^ 27(ri + r2) ^ {^i+^2Y ^ 4^ 



T^V 



112 



T1T2 



(51) 



From the result ([50| we obtain some general observations: 
The population difference goes to zero: 

• If A^ - ^1 



Ta 



j^.Then each level loses on the average as much as it gains. 

• If F — > 00. In this case the flopping between the levels is so fast that the population spends roughly half its 
time on each level, and consequently their differences are smeared out. 

On the other hand, the difference goes to the value ^ — ^ if the coupling becomes inefficient. This happens when 

• The coupling is weak, V ^ Q. 

• The decoherence rate is strong, 7 — >■ 00. 
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FIG. 1: The figures show the total population p22 + pii (solid line), p22 (dashed line), pn (dotted line) and the absolue value 
of P21 (dash-dotted line) as a function of time. In Case 1 (top left), the pumping parameters are Ai = 0, A2 = 1.0, the decay 
rates are r2 = 0, Fi = 1.0, 7 = 0.5, the detuning is a; = and the coupling V — 2.0. There is zero population in the initial 
state. In Case 2 (top right) the parameters and initial state are the same as in Case 1 but with oj = — 1.0. In Case 3 (bottom 
left) the parameters are Ai = 1.0, A2 = 1.0, F2 = 1.0, Fi = 1.0, 7 = 1.0, oj = —1.0 and V = 5.0, and the initial state has 
P22(0) — 1, with other density matrix elements equal to zero. In Case 4 (bottom right) parameters and initial state are as in 
Case 3, but with u = —5.0. All numerical values are given in terms of scaled dimensionless physical variables, as explained in 
the text. 



B. Numerical results 



The time evolution equation (47) contains 7 parameters with dimensions |i I ■ One may be eliminated by scaling 
the time rate, but the parameter space is still too large to allow a full systematic mapping. Some conclusions are, 
however, readily obtained. 



When we pose initial conditions p(0) and numerically integrate the equations (47 1 we can see the manifestations of 
the remarks above. The resulting time evolution is exemplified in Fig.l. After scaling the physical parametrs in Eq. 
(47) by a suitably selected time parameter, we have 6 remaining dimensionless variables. The numerical values used 
in the computations are these dimensionless parameters, which are given in the figures. 

Fast flopping rates V tend to make the final populations equal and the coherences zero. The decay times are 
determined by the relaxation rates, and during the time evolution the population difference tends to change sign V/j 
times; increasing the detuning, we find that w tends to decrease the amplitude of these oscillations. The effect of cu 
is, however, less than that of the coupling V. 

In Fig. 2, the functional M^ is plotted for the cases in Fig. 1. We see that in each case, the expectation value of 
fi decays smoothly and monotonically. 

Some details of these four cases chosen as illustration are given in the Appendix. As seen there, all eigenvalues 
imply damping as we expect from the physical interpretation. In all four cases, there are two purely damped modes. 
These are the ones describable as pure rate variables. We have, in fact, integrated a large number of parameter 
combinations, and we have found no anomalous behavior; all cases show smooth and regular behavior. The results 
are very similar to the ones in the cases presented here. We expect our analysis to apply to all cases. 



CONCLUSION 



The idea in our line of research has been to find expressions uniquely indicating the direction of time at each instant. 
This is particularly urgent in the field of quantum optics where the time evolution is conventionally described by 




FIG. 2: The expectation value of f2 as a function of time for the same four cases as in Fig. 1, from bottom to top Case 1 
(solid line), Case 2 (dashed line), Case 3 (dotted line) and Case 4 (dash-dotted line). In all cases, normalisation is chosen as 
(p(0)|Jl|p(0)) — 1. As the curves are very similar, they have been shifted vertically by 0.5 with respect to each other. The 
asymptotic value for the unshifted {p(i)|Sl|p(i)) when t — > oo is zero in all four cases, and the behaviour is monotonic as a 
function of time. All numerical values are given in terms of scaled dimensionless physical variables, as explained in the text. 

phenomenological master equations. In particular, the state space of interest comprises only a part of all possible 
states of the physical system. Thus an interpretation in terms of quantum concepts is natural, but also Markovian 
rate models fall under our concepts [9]. In fact, the rate equation situation is introducing an approach [18|, which 
constitutes a precursor to our method. 

In mathematics, the related formalism has been known for some time under the name of " symmetrization" |191l20j . 
The Prigogine school has introduced a similar description of irreversible time evolution |21j . but their formalism 
is based on a different physical argument and less concisely defined than our present approach. We base all our 
expressions on precisely defined mathematical concepts. When applicable our formulation leads to uniquely determined 
numerical evaluations. 

The approach is limited to linear evolution equations with time-independent generators. The linearity is natural 
if we wish time evolution in independent subsystems to remain uncoupled. Violating this condition takes us into a 
totally different category of models. To introduce generators depending on time is possible, but the ensuing theory 
then has to deal with time dependent eigenvalues and eigenelements. This is possible but leads to a theory too 
complex to give transparent physical interpretations. 

Wc feel that even with its restrictions, our method is general enough to illuminate the properties of irreversible 
system dynamics. 

Appendix A: Details of the numerical simulations 

Below we give the steady states and the eigenvalues for C relating to cases 1-4 plotted in Figs. 1 and 2. The steady 
states po agree with Eq. (50), and all real parts of eigenvalues are negative as required. The parameters are as follows: 



• Case 1: Ai = 0, A2 = 1.0, A21 = 0, Fa = 0, Fi = 1.0, 7 = 0.5, uj ^ 0, V = 2.0 
Po = (1.0625, -0.25i,0.25i,l)^ 

Ai = -0.5 - 3.9686i, A2 = -0.5 + 3.9686i, A3 = A4 = -0.5 

• Case 2: Ai = 0, A2 = 1.0, A21 = 0, F2 = 0, Fi = 1.0, 7 = 0.5, uj = -1.0, V = 2.0 
Pa = (1.3125, -0.5 - 0.25i, -0.5 + 0.25i, 1)^ 

Ai = -0.5000 - 4.0945i, A2 = -0.5000 + 4.0945i, A3 == -0.6221, A4 = -0.3779 



• Case 3: Ai = 1.0, A2 == 1.0, A21 = 0, F2 = 1.0, Fi = 1.0, 7 = 1.0, a; = 1.0, V 
Po = (1,0,0,1)^ 
Ai = -1.0 - 10.0499i, A2 = -1.0 + 10.0499i, A3 = A4 = -1.0 



5.0 



10 



• Case 4: Ai = 1.0, A2 = 1.0, A21 == 0, Ta == 1.0, Ti == 1.0, 7 = 1.0, lu = 5.0, V ^ 5.0 
Po = (1,0,0,1)^ 
Ai == -1.0 - 11.1803i, A2 = -1.0 + 11.1803i, A3 = A4 = -1.0. 
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